Kapitel IX. 


Differenzenrechnung. 
Von H. E. Timerding in Braunschweig. 


$1. Differenzen. 


Die Differenzenrechnung ist die Zwillingsschwester der Dif- 
ferentialrechnung. Die Entdecker der letzteren, Leibniz und 
Newton, haben gleichzeitig auch zuerst das Wesen und den 
methodischen Charakter der ersteren erkannt. (Vgl. Cantor, 
Gesch. d. Math. 3 (1898) S. 37, 358 und Zeuthen, Gesch. d. 
Math. im XVI. u. XVII. Jahrh.) Als eigentlicher Begründer 
der Differenzenrechnung hat Newton zu gelten, der zu ihrer 
Entwicklung in einer kleinen, 1711 erschienenen Schrift Me- 
thodus differentialis den Grund legte. (Diese Schrift ist wieder- 
abgedruckt in den Opuscula math. von Newton 1744 und steht 
in der Ausgabe von Horsley (1779) in Bd. 1, S. 519.) Als 
das erste Lehrbuch ist anzusehen Brook Taylor’s Methodus in- 
erementorum directa et inversa, London 1715, dem 1730 Stir- 
lings Methodus differentialis folgte. Eine weitere methodische 
Ausgestaltung hat die Differenzenrechnung in Eulers Institutiones 
caleuli differentialis, Petersburg 1755, erfahren. 

Die Differenzenrechnung setzt eine tabellarisch gegebene 
Funktion voraus, d.h. man nimmt an, die Werte der Funktion 
f(x) seien für eine Reihe von Werten des Argumentes x, die 
in den gleichen Abständen k aufeinander folgen, bekannt. Dann 


schreibt man: 
Afl@)=fl@+h)— f (2) 


und bezeichnet dies als die erste Differenz der Funktion f(x). 
Verfährt man mit der ersten Differenz ebenso wie vorher mit der 
Funktion selbst, so erhält man die zweite Differenz A? f(x) usw. 
Allgemein ist: 


Ar f(x) = Ar!f(a + Ah) — A" fle), 
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und man findet: 
Artnr) = ATArfle). 
Bezeichnet man der Kürze halber die Werte 
to yth mr2h... u + nh 


Zo LA Da EN E 


mit 
so wird: 
Ara) = f (8) — (1) En- + (Et) Er): 


Mit Hilfe der sukzessiven Differenzen für den Wert des 
Argumentes x, kann man den Funktionswert für x, folgender- 
maßen ausdrücken: 


fa) = fa) + (1) AFl) + (3) AE) + + Art) 
oder symbolisch: 
fa) = (1 + AY f (20). 

Die Differenz einer ganzen Funktion von v ist eine ganze 
Funktion, deren Grad um 1 niedriger ist. Die n' Differenz 
einer ganzen Funktion n*™ Grades ist somit eine Konstante, 
und zwar wird, wenn 

fl) = aa" + qaty 
angenommen wird: 
Ar f (x)= n!a,: Mh". 

Insbesondere ergibt sich: 
Ax(z—h)(2—2h)...(e—nh)—(n+1)h-2(c—h)...(«&—(n—1)h), 
außerdem die von Nicole (Hist. de V Acad. des Sciences de Paris, 
Année 1717, p. 7) gefundene Formel: 

1 1 
A- = — nh — — — ° 
æ(x +h)(æ + 2h)... (Œ+ n— 1)h) x(æ +h)... (æ+ nh) 

Die Differenzen der Exponentialfunktion sind dem Funktions- 
werte proportional. Es wird: 

Ara” = (a — 1)". a” 

Für die Kreisfunktionen findet man: 


A” sin u = (2 sin A)" -sin (u +2 (a +1), 


Ar cos u = (2 sin PDE cos (u + > (a + h)). 
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Es gelten die allgemeinen Regeln: 
Acf(x)=cAf(x) (= const), 
Alpe) + y(2)] = Ayla) + Ayla), 
Alp (2) y(æ)] = p(x) : Ayla) + yle + h) > Ag (2). 
Allgemein wird für æ + nh = %,: 
Arlp(a)- va] plr): Aula) + (1) Ap (8a —h) - Arte) 
+ (3) A?y(x,— 2h): A”-2 y (8) +++ A" g(x). y (2). 
Die Differenzen werden zu den Differentialquotienten in Be- 
ziehung gesetzt durch folgende Formel: 
A” f(x) = Ari" f(x) ++ Ar h” f™ (x) + Artiyetifatn(e), 
wo & einen Wert zwischen x und æ + nh bedeutet, ferner 
fox = a f(x) 
da“ 
und (für u > m) 


g-i [Leza 


1 
u! ht P (Are); 


z=a 


en a A A 
gesetzt ist, Umgekehrt wird: 
Kaf (4) = BZA” f (2) +: + BEA" (a) + Buttmntigatn (E), 
Hierin ist & wieder ein Wert zwischen æ und æ + nh, und es ist: 


air la rag 1 


Die Koeffizienten Ah, und B}, kann man auch durch die Reihen- 
entwicklungen festlegen: 


(è — 1)" = Ah” + Amtipmtihe.., 
Dog (1 +)" = Bam + Brttimtt L..... 


Daraus folgt die nachstehende symbolische Schreibweise der 
obigen Gleichungen: 


are ine, Ge". 


Pascal, Repertorium. I. 2. Aufl. 33 
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Die Gleichungen sind so zu verstehen, daß der symbolische 
Faktor von f (æ) auf der rechten Seite nach Potenzen von D, 
resp. A entwickelt und hierauf 


Di, fa) = En, Au. f(a) = A" f(@) 


gesetzt werden muß. 

Mit Hilfe dieser symbolischen Schreibweise kann man leicht 
eine Formel geben, die es gestattet, von den für ein Intervall A 
gebildeten Differenzen A“f(x) zu den für ein anderes Intervall 
h, gebildeten Differenzen A”f(x) überzugehen. Diese Formel 
lautet: 


Ant) = [1 + A)" — 17" f(a). 


$ 2. Interpolation. 


Die Aufgabe der rationalen Interpolation verlangt zunächst, 
eine rationale, insbesondere eine ganze Funktion f(x) so zu be- 
stimmen, daß sie für gewisse Werte des Argumentes 


Ua e Da 


vorgegebene Werte 
Yo Y1 Ya - -+ Yn 


annimmt. Es sind zwei verschiedene Lösungsformen dieser Auf- 
gabe zu unterscheiden, die Lagrangesche und die Newtonsche 
Die Lagrangesche Formel erhält man, wenn man 


p(z) = (z — t9) (@ — 2) . . . (€ — 2p)» 
d 
ee... en e) 
setzt, in der Gestalt: 
uhr, 
n) fi PR) E— Lu 


u=0 


Die Newtonsche Formel lautet in der Fassung, die ihr 
Gauß gegeben hat: 


f (£) = Y + (£ — z0) [£0 21] + @— zo) (® — 21) [29 £1 to] 


+(x — 2) (£ — 2) (@ — 2) [t02 % slt: 
Hierbei ist: 
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=N Ya Yy—Y% 
Mal le et N ee 
[a 2] — [2o 2] [2; xs] — [2 2] 
[22] = a E, lera a A A 
[X &; &] — [2 2 2] 
[ro ti 2%] = = a PAS» 


Die Newtonsche Formel hat den Vorzug, daß man bei 
Hinzunahme neuer gegebener Größen nur neue Glieder hinzu- 
zufügen braucht, ohne die vorhergehenden ändern zu müssen. 

Die weitere Aufgabe ist nun die, von einer bestimmten 
Funktion F(x), deren Werte für eine Reihe von Argumenten- 
werten gegeben sind, die Werte für alle Argumentenwerte eines 
gewissen Bereiches näherungsweise zu ermitteln. Die vorstehend 
angegebenen ganzen Funktionen f(x) sind dann als Näherungs- 
formeln für die Funktion F(x) anzusehen, und es wird somit: 


F(z) =f(2)+ R, 
wo R ein Korrektionsglied bedeutet. Nach Cauchy, ©. R. 11 
(1840) 787, (Euvres (1) 5, 422, hat man zu setzen: 
—%)@— 2) . . (C—&n) sing 
EEE er), 
wenn & einen Wert zwischen «, und x, bedeutet. 

Für den besonders wichtigen Fall, daß die Werte der Argu- 
mente, für welche die Funktionswerte gegeben sind, in gleichen 
Intervallen % aufeinander folgen, also von der Form sind: 

Zos To F h, £o + 2h o.. ao + nh, 
geht die Newtonsche Formel über in: 


Oe ee A Are 2 @— en £o —h) pa "22 


en re Afa) 
% n! ET RER 


Substituiert man: 
æ = To + ht 
und setzt: 
(= + hrt) = g(t), Fl) = PO) = v, 
so wird einfacher: 
pÀ =n FAm aia TR E u Pe er 
33* 
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Diese Funktion nimmt für 
reo B, 


die verlangten Werte an, und so wird die Newtonsche Formel 
gewöhnlich gegeben. 

Durch größere Symmetrie ausgezeichnet ist die Formel, die 
dazu dient, die Werte der Funktion in der Umgebung einer be- 
stimmten Stelle zu bestimmen und die Stirlingsche Formel 
genannt wird, obwohl sie in einem von Newton gegebenen 
Satze (Methodus differentialis, Prop. IIT) bereits enthalten ist. 
Wir nehmen auch die Werte des Argumentes t=— 1, — 2... 
hinzu und stellen die folgende Tabelle auf: 


Arg. | Fkt. Differenzen 
t y I Ti II IV N YI 
— 3 y | 1 
aA a E AE TARI 
Dea ag | A| a A n au 
0 Yo Ay, A°y, A? Fr Ay, A5 7 A Ys 
1| 9 Ay | ya | 9 
Ay, 2y | A'Y i 
: Ya Ay, A®y, 
Ys | 


Setzen wir dann noch in leicht verständlicher Symbolik: 
An + Aw) AI, IA’ + A’) = A... 


so wird die Stirlingsche Formel, indem t einen Wert zwischen 
— 4 und + 4 bedeutet: 


g(r) =% + TAJ, +3, Ad + a leeren) Ay, 
ER u z DNE +23 + ar New 2A, u 


Die wichtigste Literatur über Interpolation sei im folgenden 
zusammengestellt: Lagrange, Sur les interpolations, Œuvres 7, 
535; Sur une methode partic. d’approximation, (Euvres 5, 517; 
Mém. sur la methode d’interpolation, CEuvres 5, 663; Sur une 
nowvelle espèce de calcul, (Euvres 3, 441. Gauß, Theoria inter- 
polationis, Werke 3, 265. Encke, Über Interpolation, Berl. 
astron. Jahrbuch 1830; Über mechanische Quadratur, ebenda 
1837, 1862; Ges. Abh. 1, 1, 21, 61. Cauchy, Mém. sur Vin- 
terpolation, J. de math. 2 (1837) 193. F. Neumann, Astr. Nachr. 
15 (1838) 313. Villarceau, Add. à la Connaiss. des Temps 
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1852, p. 129. Genoechi, Ann. di mat. (1) 6 (1864). Tsche- 
byscheff, Petersburg. Mém. (T) 1 (1859); Œwvres 1, p. 387, 
und ebenda p. 539. Hermite, C. R. 48 (1859) 62; J. f. Math. 
84 (1878) 70. Merrifield, Report Brit. Assoc. 1880. Gram, 
J. f. Math. 94 (1883) 41. Méray, Ann. Ec. norm. (3) 1 (1884) 
165. Harzer, Astr. Nachr. 115 (1886) 337. Radau, Études 
sur les formules d’interpolation, Paris 1891. Netto, Math. Ann. 
42 (1893) 453. Pincherle, Bologna Mem. (5) 3 (1893) 293. 
Tisserand, Mécanique céleste (1896), Vol. 4, p. 15 (Leverriers 
Interpolation durch periodische Reihen). Bruns, Astr. Nachr. 
146 (1898) 161. Runge, Ztschr. f. Math. 48 (1903) 443. 
De la Vallée Poussin, Brus. Acad. roy. Bull. (1908) 319. 
Schließlich sei das kürzlich erschienene Lehrbuch von Thiele, 
Interpolationsrechnung (1910) genannt. 


§ 3. Summen. 


Zu der Differenz invers ist die Summe. Wenn die Funktion 
f(x) gegeben ist, so sucht man die Funktion g(x), welche der 
Funktionalgleichung genügt: 


g(a + h) — g (2) = f (2), 


deren Differenz also die gegebene Funktion ist. Durch die vor- 
stehende Gleichung wird die Funktion: 


g (2) = Url) 


für eine Reihe von Werten x, die in dem Intervall } aufein- 
anderfolgen, bis auf eine willkürliche Konstante, die Summations- 
konstante, bestimmt. Für die Summen gelten analoge Glei- 
chungen wie für die Differenzen. Es ist: 


Dofl@)= CHE) (C= const), 
Dlol) + ya] = Zol) + Ey), 
außerdem gilt die Formel der partiellen Summation: 
pl): Ayla) = Pla) ya) — Dye + h) - Ayla). 


Ferner wird, bis auf die noch hinzutretende Summations- 
konstante: 


De-i)... @-m—DN-„np2@—h) ... (&— nh). 


Allgemein ist die Summe einer ganzen Funktion n'" Grades 
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f„(x) eine ganze Funktion (n + 1)® Grades f,,.(x). Ist A,2” 
a+ 
das höchste Glied in f (x), so beginnt f,,,(x) mit A, 
Endlich gelten die Formeln: 


J 1 T 1 1 
Pa sæ +h)... (@nh) nhaac-+h)...etm—ın’ 


1 
> ANZ ar 
1 
co: & si > 
f s(u-5) in (3) 
X'sinu=— —— 4, Fou = — "7; 


N, EN N 
2 Ga k ben 
sın 3 2 sin 2 


Onti 


Von den unbestimmten Summen, die noch eine willkürliche 
Konstante enthalten, geht man zu den bestimmten Summen über, 
indem man die Konstante dadurch festlegt, daß für einen be- 
stimmten Wert des Argumentes die Summe verschwinden soll. 
Dieser Wert x, heißt die untere Grenze der bestimmten Summe, 
und man schreibt diese in der Form: 


z 
PAGE 
Aus der vorangestellten Definition des Begriffes der Summen 
ergibt sich dann sofort, daß: 
Zo +nh 
2i fC) = f (€) + fC tR) t: flo + an DA 


ist. Nimmt man z. B.: 
f(&) 8, 


so erhält man: 


mi (x — h) i R i 
> t y (Summenformel der arithmetischen Reihe); 
0 


nimmt man: 
f(z) =", 


so wird: 


x x 
wi (Summenformel der geometrischen Reihe). 
0 a1 


Für f(x) = x° ergibt sich die Summe: 


Deal r Tela 


Der wesentliche Inhalt dieser Formel war nach M. Cantor 
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(Gesch. d. Math. 1, 724) schon im 11. Jahrhundert den Arabern 
bekannt. 
Man findet ferner: 


sin (« z) sin h 
t EELA De 
nn > cos u = RASTE ANTY 
h 2 sin ar A 
2 


nR 
c08 — — COS |U — — 


u 
5, sinn La un 2 
h ` h 


2 sin a 


Wir wollen nun die Summen 


= x” 
GE DET, 


betrachten, die nach J. L. Raabe-als Bernoullische Funktionen 
bezeichnet werden (Die Jacob-Bernoullische Funktion. Zürich 
1848). Nach dem, was wir oben über die Summe einer ganzen 
Funktion nt Grades gesagt mak können wir schreiben: 


-1 
h Pay) = Er tAk E Ah? ot + AME. 
Da: 


Pn+3(2) = nu de, 


MEPA sich, wenn man von P,+ı(®) zu @,;s(®) übergeht, die 
Konstanten A,... A, nicht, es tritt nur eine neue Konstante 
A,;ı hinzu. Für = h wird Pn+1(8) = 0; somit ergeben sich 
zur Bestimmung der Koeffizienten die Gleichungen: 


1 


1 1 1 
y un m 


4=0,4 0 usw. Jedes A mit ungeradem Index, 


So. findet man: A, = — 4=-0,4,=- 


ausgenommen A,, verschwindet. Die Funktionen 9,,(x) mit 
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geradem Index enden also mit dem Glied, das gs? enthält, 


Ferner wird: 
Por (h — 2) = pa, (2); 
dagegen: } 
h 
Por (AR — 2) = — Paz+1 (2), woraus Par+ı (z) AP 


Die A mit geradem Index, die außer A, allein übrig bleiben, 
sind abwechselnd positiv und negativ. Es wird gesetzt: 


A=-B, ud 4=- (k= 2,3...) 


also 


B,=4; B=-B=-B,=::-0. 
Für B,, erhält man dann die Reihe: 


k)! 
(2k) ER E ak a i 


T (- 12 amt 


Diese Entwicklung hat Euler gegeben in den Comment. 
Acad. Petropol. 12 2 p- 73. Man vergleiche auch Schlö- 
milch, Ztschr. f. Math. 1 (1856) 200, Sonin, J. f. Math. 116 
(1896) 138. Die Zahlen B,, heißen Bernoullische Zahlen. 
Sie können durch die symbolische Gleichung 


(B+1"— B=n (n=1,3,3,...) 


definiert werden, in der nach der Ausrechnug B* durch B, zu er- 
setzen ist. Die ersten unter ihnen sind: 


1 1 1 1 5 
Bug Bug Ben Ben zo Ba ano 
691 7 3617 43 867 174611 
Bas = — z0 Pu = g> Pi = go: Bis -aa ı Pom -ao 


Nach von Staudt (J. f. Math. 21 (1840) 372) ist 


1 
Bu=n- 7: 


wo n eine ganze Zahl und die Summation nach allen Prim- 
zahlen l zu nehmen ist, für die Z—1 ein Teiler von 2% ist. So wird 


Die ersten 62 von Null verschiedenen Bernoullischen Zahlen 
hat Adams, J. f. Math. 85 (1878) 269, die ersten 90 Ssere- 
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brennikow, Petersb. Abh. (8) 16, No. 10 (1905), berechnet. Lite- 
ratur über diese Zahlen bei Ely, Amer. J. of Math. 5 (1882) 228. 

Aus der Definitionsgleichung der Summe entsteht eine 
Differenzengleichung, wenn die Differenz der gesuchten Funktion 
y = p(x) nicht unmittelbar als eine bekannte Funktion von x, 
sondern als Funktion von x und y zusammen gegeben ist, ins- 
besondere als eine für y lineare Funktion: 


Ay Fr P,y ar Q3 
wo P,, Q, bekannte Funktionen von æ bezeichnen. Schreibt man: 
Y=Ys DY = Yey Ya 


indem man das Intervall der Argumente h = 1 annimmt, so 
wird die vorige Gleichung von der Form: 


Yzy FY A,Y, 225 Qo- 
Allgemein heißt jede Gleichung von der Form: 


Yrm A A ian an Anz tr ET ar Am -dg = Q, 
eine lineare Differenzengleichung m‘ Ordnung, indem man als 
Differenzengleichung schlechthin jede Gleichung zwischen den 
Werten des Argumentes æ und den Funktionswerten Yp, Ypi 
für eine von x anfangende Reihe äquidistanter Argumentenwerte 
bezeichnet. 

Näheres über Differenzengleichungen siehe in Kap. X. 


$ 4. Zusammenhang zwischen Summen und Integralen. 
Mechanische Quadratur. 


Wie wir oben von dem Zusammenhange zwischen Differenzen 
und Differentialquotienten gesprochen haben, so kann man auch 
eine Beziehung zwischen Summen und Integralen suchen. 

Auf der einen Seite kann man die Summe durch ein In- 
tegral mit Zuhilfenahme sukzessiver Differentialquotienten aus- 
drücken. Dieser Ausdruck wird durch die Eulersche Formel 
geliefert: 


h ar (2) s Stæ da + A,h[f(8) — fla) + Air (0) — f’(a)] 
a e T O a E 


A,, A ... bedeuten hierin die früher so bezeichneten Zahlen, 


und es ist h = bea Das Restglied R,,, hat Poisson, Par. 
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mem. 6, 580 (1826) hinzugefügt; es wird nach Jacobi (J. f. 
Math. 12 (1834) 263, Werke 6, 64) in folgender Form dargestellt: 


h 
b 
Brp = — fon du D fOe +R — u), 
0 z 


wenn p„}1(%) wie oben die Bernoullische Funktion bedeutet, 
und es wird 


Bus Aaa E T E a E D): 


Die Formel ohne Restglied wurde von Euler gegeben, Comment. 
Acad. Petrop. 6 (1738) 68 und ebenda 8 (1741) 147. Man 
sehe auch Malmstén, Acta math. 5 (1884) 1, Sonin, Ann. ee. 
norm. (3) 6 (1889) 257, C. R. 108 (1889) 725, Kronecker, 
Vorl. üb. die Theorie d. Integrale (1894), 130—156. 
Als erstes Beispiel nehmen wir: 

f(&) = =", b=x, a=0, 
so wird: 

x am+1 Di 


Nimmt man zweitens: 
f(x) =e bel, E A 
so wird: 1 
e—ı 


- + tt Arte. 


Nimmt man: 
f(z) =log&, =m, a=; k=l, 


und addiert zu beiden Seiten der Gleichung log m, so erhält 
man die Stirlingsche Formel (Stirling, Methodus differentialis, 
1730, p. 135): 


log (m!) = log 2m + (m+ 5) log m — m + 4, +e 


2k — 2)! 2k)! 
+ Ar IF 0 darsa T 
(0< 0zi 


Auf der anderen Seite kann man nun auch näherungsweise ein 
Integral durch Summen darzustellen suchen. Diese angenäherte 
Auswertung heißt mechanische Quadratur. Die einfachsten Formeln, 
die dies leisten, sind die Trapezformel und die Simpsonsche Regel. 

Man denkt sich das Intervall von a bis b in n gleiche 
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Teile h zerlegt, so daß b — a = nh, und bezeichnet die Funk- 
tionswerte an den Grenzstellen x, = a + uh der Teile mit: 


Yos Yir Ya - + + Yni 
die Fumktionswerte in der Mitte der Teile durch: 


Y» Yy ee 


Dann wird die Trapezformel: 


fioa =E (nt 2n 2m H H 2ni FR 


mit: 


R, nr), a<t<b. 
Die Simpsonsche Regel lautet: 


b 
rn 2y,+4yy +49, E 


Simpson, Math. diss., London 1743, jedoch bereits bei Gregory, 
Exerc. geom., 1668 (vgl. Heinrich, Bibl. math. (3) 1, 90 (1900)) 
und Toricelli, Quadrat. parabolae, 1644 (vgl. Aubry, J. de 
math. élem. (5) 21 (1897)). 

Die erste Formel, welche die angenäherte Quadratur in 
konsequenten Zusammenhang mit der Theorie der Interpolation 
brachte, war die Cotessche Formel. In ihr wird das zu berech- 
nende Integral ersetzt durch das Integral einer ganzen Funktion 

ten Grades, welche für x£), %4 -.. 8, dieselben Werte Yo, Yı - -- Yn 
annimmt wie die zu integrierende Funktion. Setzt man: 


g(x) = (£ — 20) (£ — 2) -© (£ — Tp) 


und schreibt: 


b 
0 1 n 
Siou- (b — a) [Hy -F Hy SEE Hyn], 


indem man sich f(x) durch die Lagrangesche Interpolations- 
formel bestimmt denkt, so hat man zu setzen: 
b 
4 (æ) dæ 
H=-|- a: 
n ee (b — a) p (£) 
Substituiert man hierin: 


www.rcin.org.pl 


524 Kapitel IX. Differenzenrechnung. 


x=a-+hz, =athu, 


so kann man das Integral von der besonderen Wahl der Grenzen 
unabhängig machen und schreiben: 


a mar B-etDe en. om, 


e(w— I). IE er (a —n) 


Hierbei wird: # De S je d 
n n u=0 n 


Wir stellen in der folgenden Tabelle die von Cotes (Har- 
monia mensurarum, 1722) berechneten Werte für n = 1 bis n = 10 
zusammen: 


u n-u| u 
H=H 
nn 0 1 2 3 4 5 
n—1 1 | i 
2 2 
n= 2 4 4 1 
6 6 6 
n = 8 1 3 2 TA 
8 8 8 8 
EN a 7 
90 90 90 90 90 
| A 5 0 50, BR 
288 288 288 288 280 288 
> 216 21 272 21 216 
840 840 840 840 840 840 
San 751 3577 1323 2989 2989 1323 
17280 | 17280 17280 | 17280 17280 | 17280 
a We... 5888 928 10496 | 4540 | 10496 
28350 | 28350 | 28350 | 28350 28350 | 28350 
ng | 2857 15741 1080 | 19344 56778 | 5778 
89600 | 89600 89600 | 89600 89600 | 89600 
n — 10 | „16067 | 106300 48525 | 272400 260550 | 427368 
598752 | 598752 598752 | 598752 598752 | 598752. 


www.rcin.org.pl 


u 


§ 4. Mechanische Quadratur. 525 


Gauß hat den Gedanken verfolgt, die Werte des Argu- 
mentes £y, Ty ... Lpy welche die Grundlage der Interpolation 
bilden, nicht durch Teilung des ganzen Intervalles in gleiche 
Teile zu gewinnen, sondern sie so zu wählen, daß die Genauigkeit 
der Formel eine möglichst große wird. Dies sucht er zu erreichen, 
indem er die Forderung aufstellt, daß die Formel nicht bloß 
exakt sein soll für die Funktion n‘" Grades, welche für £g, &, ... & 
die vorgeschriebenen Werte annimmt, sondern für jede ganze 
Funktion bis zum (2n + 1)‘ Grade, welche dieser Bedingung 
genügt. Das Schwergewicht liegt somit in der Festlegung der 
Funktion: 


9a@)=-@-n)e— 3) (Œ — 2); 


welche, gleich Null gesetzt, die Argumentenwerte &,...%, 
liefert Das Resultat ist, daß der aufgestellten Forderung ge- 
nügt wird, wenn: 


pla) = ale att (e — b)+] 
angenommen wird. In der so entstehenden Gleichung p(x) = 0 


werden die von Fall zu Fall wechselnden Grenzen a, b zu den 
festen Werten 0,1, wenn 


| 2=a+(b—a)t 
substituiert wird. Dann wird die Gleichung: 
n+1 
Taina- y= o, 


und ihre Wurzeln können ein für allemal berechnet werden. 

‘ Diese Berechnung hat Gauß bis auf 16 Dezimalen durchgeführt. 
Die folgende Tabelle gibt die Wurzeln für die ersten 7 Grade 
auf 6 Dezimalen: 


Grad Wurzeln t 


n+i | 

T ty t ts ts t, A ts r 
E = 

1 0,5 | | 

2 0,211325 | 0,789675 | | 

3 0,112702 | 0,5 0,887298 | 

4 0,069432 | 0,330009 | 0,669991 | 0,930568 

5 0,046910 | 0,230765 | 0,5 0,769235 | 0,953090 

6 — | 0,033765 | 0,169395 0, 380690 | 0,619310 | 0,830605 | 0,966235 

7 0,025446 | 0,129234 | 0,297077 0,5 0,702923 | 0,870766 | 0,974554 
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Es sind nun weiter zu berechnen die Integrale: 


ad tern 4) 
OETAN w tr e Eut 


Die Werte derselben für die ersten 7 Grade sind in der fol- 
genden Tabelle zusammengestellt: 


Grad 
1 | | x x x 
nE E, E, | E, E, E, E; E, 
1 1 
1 
A | 
2 Bro) 
š 5 TENT.: 
18 18. 31.138 | | | 
4 | 0,173927 |0,326073| 0,326073 | 0,173927 | | 
5 | 0,118463 |0,239314 | 0,284444 | 0,239314 | 0,118463 
6 | 0,085662 |0,180381 | 0,233957 | 0,233957 | 0,180381 |0,085662 
7 | 0,064742 | 0,139853 0.190915 | 0,208980 | 0,190915 0,139853 0,064742 


; Die Gaußsche Formel schreibt sich dann in der einfachen 
Gestalt: 


Ji@dz=(b0—a)[E f(a + 0—0) +-+ Efla + 0—0). 


Außer der Arbeit von Gauß (Comment. Gotting. 3 (1814/15), 
Selbstanzeige Gött. gel. Anz. 1814, Werke 3, 163, 202) ver- 
gleiche man u. a.: Jacobi, J. f. Math. 1 (1826) 301, Werke 
4, 3. Christoffel, J. f. Math. 55 (1858) 61. Tscheby- 
scheff, Journ. de Math. (2) 19 (1874) 19. Mansion, Mathesis 
1 (1881), Supplem., p. 1; Bruxelles Ann. Soc. Scient. 5 B 
(1881) 231; Bruxelles Bullet. (3) 11 (1886) 293; Comptes 
rendus 102 (1886) 412 (Bestimmung des Restes in der Gauß- 
schen Formel). Stieltjes, Ann. Éc. norm. (3) 1 (1884) 409; 
Comptes rendus 99 (1884) 850. Markoff, Math. Ann. 25 
(1885) 417. Man vgl. auch den Bericht von Burkhardt, 
Math.-Ver. 10 2, § 88. 


Lehrbücher der Differenzenrechnung: 


S. F. Lacroix, Traité des différences, Paris 1800, 3. Bd. 
des Traité du calcul differentiel etc. J. F. W. Herschel, Collec- 
tion of examples of the applications of the calculus of finite diffe- 
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rences, Cambridge 1820. O. Schlömilch, Theorie der Diffe- 
renzen und Summen, Halle 1848. G. Boole, A Treatise on 
the calculus of finite differences, Cambridge 1860, deutsch von 
Schnuse, Braunschweig 1867 u. d. T.: Die Grundlehren der 
endlichen Differenzen- und Summenrechnung. A. Markoff, 
Differenzenrechnung, deutsch von Friesendorff und Prümm, 
Leipzig 1896. E. Pascal, Calcolo delle differenze finite, Milano 
1897. H. Bruns, Grundlinien des wissenschaft. Rechnens, Leipzig 
1903. O. Biermann, Vorlesungen über mathematische Näherungs- 
methoden, Braunschweig 1905. D. Seliwanoff, Differenzen- 
rechnung, Teubners Sammlung Bd. 13, Leipzig 1904. 


Verbesserungen. 
S, 147, die letzten vier Zeilen sind zu streichen. 
„ 213, Zeile 22 v. o. lies: |z, «+ ße | statt |z, az +£]. 
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